function orthbasis(t)

if nargin<1
    t=1;
end
ms{1}=@(s) 1/sqrt(t)+s-s;

for kk=2:100
    ms{kk}=@(s) sqrt(2/t)*cos((kk-1)*pi*s/t);
end

fname='orthbase';
save(fname,'ms');

%
% for ii=1:10
%     for jj=1:10
%         disp(quadl(@(s)ms{ii}(s).*ms{jj}(s),0,1));
%     end
% end